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Abstract 


A probabilistic risk assessment procedure is followed to estimate the probability of protection at facilities 
from vehicle-borne explosions. Using truncated normal distributions for both the TNT equivalencies of 
the charge and the standoff distances, a probability for protection is calculated for a single degree of 
freedom (SDOF) blast protection design curve. Numerical examples are given for two representative 
design curves that give probabilities of protection of 80% and 91%. 
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1. Introduction 


Different scenarios can be envisioned that would be relevant to the design of terrorist mitigation measures 
at a facility. For explosive threats, one could assume a scenario in which the standoff distance is the 
minimum possible distance to the intended target, produces collateral damage to nearby facilities and that 
the explosive source is very large. For example, one such scenario might be 40,000 lb of an efficiently 
exploded ammonium nitrate/fuel oil mixture contained within an 18-wheeler. A related scenario might be 
where the 18-wheeler is located as close as physically possible to the intended target or where the size of 
the explosive source is not as large. If a smaller vehicle is envisioned to transport the source, 40,000 lb 
would be an overestimate of the explosive payload. In addition, the TNT-equivalent mass of the 
explosive can vary depending on the physical and chemical details of the explosive mixture. 


The uncertainties in explosive strength and standoff distance suggest that a statistical assessment of the 
threat parameters could play a role in the mitigation design. This type of approach has been proffered 
recently by Stewart, Netherton, and Rosowsky (2006). In this reference they argue that there is so much 
uncertainty associated with terrorism that an attempt should be made to quantify the uncertainties, and 
that these results should then be used in established probabilistic risk assessment procedures to 
systematically assess the viability and relative benefits of different mitigation measures. 


Netherton and Stewart (2009) considered the variability in explosive blast loading caused by a 
multiplicity of factors with respect to facade glazing, and used Monte Carlo simulation to calculate 
probabilities of glazing damage and safety hazards conditional on scenario-driven parameters. 


The present paper also provides a statistical assessment of an explosive-borne threat to window structures 
in facilities. However, we adopt an exclusively analytic approach based on published scaling relations for 
the key explosive parameters, impulse and overpressure, as well as the design specifications of the 
window. Although it does not take into consideration every physical parameter that affects explosive 
impact (most notably the presence of intervening structures), this method obviates the need for computer 
simulation and yields a simple if approximate method to assess explosive risk. 


Specifically, we estimate the probability of protection one might expect against likely blast overpressures 
and impulses. This is based on both the probability function for explosive strength and standoff distance 
and the design curve used for constructing components of a curtain wall facade, a popular form of 
construction for modern office facilities. In this paper, the design curve is derived from a simple single 
degree of freedom (SDOF) component model. 


In addition, the method of assessing the probability of protection afforded by a specific building 
component as described herein can be generalized and applied to other security scenarios for which 
scaling relations exist for relevant physical parameters. This offers a potentially powerful tool in 
quantitatively assessing risk for a variety of scenarios. 


Section 2 describes the probability distributions to be assumed for the TNT equivalence strength of the 
explosive source, and for the standoff distance between the explosive source and a potential target. 


Section 3 discusses the blast parameters (overpressure and impulse) that determine the damage 
experienced by a structure. 


Section 4 summarizes scaling laws that relate the blast overpressure and impulse to the explosive strength 
and standoff distance. 
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Section 5 then discusses a blast protection design curve that is based on a SDOF model and on data 
provided by a consulting firm for a specific window design. 


Section 6 derives a specific blast mitigation probability of protection using the probability function and 
the design curve. 


Section 7 discusses the specific results and generalizes the approach to other operational risk problems. 


2. Probability distributions for explosive source strength and for standoff distance 


Based on past terrorist explosive attacks on structures, we shall assume that the explosive source is an 
ammonium nitrate/fuel 011 (ANFO) mixture, and that this explosive mixture is delivered to the near 
vicinity of a facility by a vehicle. Uncertainties exist as to both the TNT equivalency of the ANFO, and 
to how close the vehicle will be to a facility when the mixture is detonated. 


2a. Uncertainty in TNT equivalency factor 


The construction and the efficacy of the resulting bomb have been widely discussed in the open literature. 
Two readily available on-line summary articles are: 


Explosives - ANFO *Ammonium nitrate-fuel oil), at 
http://www.globalsecurity.org/military/system/munition/explosives-anfo.htm 


ANFO, at http://en.wikipedia.org/wiki/ANFO 


Differences in the preparation of ANFO explosive as well as the physical make-up of its components can 
introduce variability in TNT equivalency. These can be manifest as variability in the amount of absorbed 
water in the mixture, a lack of uniformity in mixing, differences in specific gravity, and the particular 
oxidizer and type of fuel oil used. 


Because of the inherent variability in one or more of these factors, there is uncertainty as to the TNT 
equivalent mass that should be assumed for a terrorist attack. To provide a framework for assessing the 
blast parameters, we shall assume in the following that the probability F,,(m)dm that a terrorist explosive 
attack will involve an explosive of TNT equivalent mass m in the interval (m, m+dm) has the shape of a 
normal probability distribution function with a mean m, and a dispersion 6m, We shall also assume that 
the truncated normal probability distribution is only nonzero for m larger than some minimum mypin and 
for m less than some maximum Max. 


Specifically, 
Cn exp|[- (m-m,)°/(5m)]Jdm for Min <M < Max 
Fn(m)dm = 
0 otherwise [la] 
where the normalization constant is 


C= (2/r0'78m) [Erf((mmax-M,)/5m) + Erf((mo-mnin)/SmM) [1b] 


From the foregoing discussion, we shall use as a numerical example 
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my = 32,000 Ibs TNT 

Mmin = 12,000 lbs TNT 

Mmax = 64,000 Ibs TNT 

Sm = 20,000 lbs TNT [2] 


The mean value, 32,000 lbs, corresponds to 40,000 lb of ANFO with a TNT equivalency factor of 0.8, 
whereas the minimum and maximum values, 12,000 and 64,000 Ibs, correspond to 40,000 lb of ANFO 
with TNT equivalency factors of 0.3 and 1.6, respectively. The dispersion 6m = 20,000 lbs has been set 
equal to My — Mmin to describe a likely spread in values comparable to the difference between the nominal 
0.8 equivalency factor and the minimum equivalency factor of 0.3. This spread has been chosen rather 
than a spread equal to the difference between the maximum equivalency factor of 1.6 and the nominal 
factor of 0.8, to weight more heavily the multitude of factors that can decrease the TNT equivalency. 


The truncated normal probability distribution F,,(m) is displayed in Figure 1 for the numerical parameters 
of eq. [2]. 


Fn(m) 


0.00008 
0.000025 
0.00002 
0.000015 
0.00001 
510°° 
10000 20000 30000 40000 50000 60000 70000 
TNT equivalent mass (lbs) 


Figure 1. Probability function F,,(m) for TNT equivalent mass m for the representative parameters of eq. [2]. 


2b. Uncertainty in Standoff Distance 


Vehicle-borne explosive attacks represent a significant security concern and collateral damage 
experienced as a result of attacks against nearby facilities is a concern as well. This adds to the 
uncertainty associated with urban terrorist threats and thereby affects mitigation strategies under 
consideration. Strategically placed bollards will prevent a vehicle from ramming a facility as well as set a 
minimum vehicle access distance. We assume the closest distance of approach to a facility and still be 
within the iconic “orbit” of its intended target is 275 feet. Therefore, the standoff distance for a vehicle- 
borne explosive source should equal or exceed a 275 foot minimum. 


However, a number of vehicle positions that border a target facility are equally likely from a terrorist’s 
perspective since any point along this border would inflict approximately the same damage to the target 
facility. However, differences in this distance could have a profound effect on collateral damage to 
neighboring facilities. The range of likely detonation positions relative to potential explosive damage to a 
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neighboring facility is assumed to be principally dictated by the length of a street bordering the target 
building 


Without further details of the terrorists’ plans, we again resort to the normal probability distribution 
function as the one that best reflects the inherent uncertainty in key physical parameters. It should be 
noted that although a truncated normal distribution is most appropriate in this instance, the approach 
specified herein does not depend on the specific probability distribution. Since F,,(m) and F,(r) are 
assumed to be independent, the joint probability and resulting contour plot could be calculated in similar 
fashion using alternative distribution functions. However, the dispersion associated with the chosen 
probability distributions which characterizes the uncertainty in potential threat scenarios is considered 
most important in calculating the probability of protection afforded by a given window design. 


Specifically, we assume that the probability F,(r)dr that the standoff distance is between r and r+dr is also 
given by a truncated normal distribution: 


C, exp[- (r-to) (Sr) dr for Min < T < Max 
F,(r)dr = 
0 otherwise [3a] 


where the normalization constant is 
C.= (2/0 Sr) [Erf((tmax-To)/Sr) + Erf((to-tmin)/Sr) [3b] 


Here r, denotes the value of the standoff distance at which the probability is largest: for example, when 
the detailed geometry of a target facility placement becomes available, this could be determined from the 
location of a vehicle that would place it as close as possible to the target. The dispersion dr measures the 
spread in likely values for the standoff distance. And as discussed earlier, rmin denotes the closest possible 
distance of approach to a neighboring facility and r,,, denotes the maximum standoff distance from the 
neighboring facility for which an explosive-laden vehicle could do substantial damage to the target 
building. 


To illustrate the approach, in the following we shall present a numerical example for 
Imin= 275 feet 
Imax= 525 feet 
To = 0.5 (tmintTmax) = 400 feet) 
Sr = Lo-Tmin= 125 feet [4] 


As more detailed information becomes available, these numbers can be modified. 


The probability function F,(r) is displayed in Figure 2 for the numerical parameters of eq. [4]. 
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Figure 2. Probability function F,(r) for standoff distance r from a facility in proximity to the intended target for the 
representative parameters of eq. [4]. 


It is also of interest to look at the joint probability distribution function F,,,(m,r)dmdr that gives the 
probability that the TNT equivalent mass m is in the range (m, m+dm) and the standoff distance is in the 
range (r,rt+dr). Since we have assumed that the probability distribution functions of equations [1] and [2] 
are independent, we can write simply 


Fmr(m,r) = Fn(m) F,(r) [5] 


The joint probability function F,,,(m,r) is displayed in Figure 3 for the numerical parameters of eqs. [2] 
and [4]. As expected, the joint probability distribution function displays a single maximum at (mo, fo) 
with a spread in the TNT-equivalent mass m determined by 5m and a spread in the standoff distance r 
determined by 6r. 


Probability 


Figure 3 Joint probability function F,,,(m,r) for the numerical parameters in eqs. [2] and [4]. In the 3D image of 
Figure 3a, the x-axis displays the TNT equivalent mass in lbs, and the y-axis displays the standoff distance in feet. 
In the contour plot of Figure 3b, the TNT equivalent mass is shown on the vertical axis and the standoff distance is 
shown on the horizontal axis. 
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The contour diagram of the probability function is shown in Figure 4. 


TNT equivalent 
mass m (lbs) 


300 350 400 ASO 500 


Standoff distance r (feet) 


Figure 4. Contour plot of joint probability function F,,,(m,r) for the numerical parameters in eqs. [2] and [4]. The 
TNT equivalent mass is shown on the vertical axis and the standoff distance is shown on the horizontal axis. 


Finally, we note that a joint probability distribution could be calculated in this way for an unlimited 
number of physical parameters with an associated dispersion specified for each, assuming such 
parameters were independent. In this report, standoff distance and TNT equivalent mass were used since 
these directly affect overpressure and impulse as indicated in referenced scaling relations. 


3. Structural damage from an explosive blast 


In the preceding section we have discussed the probability distribution function Fy,(m,r)dmdr for the 
TNT-equivalent mass m and the standoff distance r. The equivalent mass m and standoff distance r are 
not the explosive blast variables that are most directly related to structural damage. Rather, these 
variables are the blast overpressure p and the blast impulse I [See, e.g., M. Held, 1983]. In Section 4 we 
shall discuss how p=p(m,r) and i=i(m,r) are related to m and r. In this section, we shall briefly review why 
it is that blast overpressure p and blast impulse i are the quantities of most direct interest to determining 
structural damage. 


The interaction of a blast wave with a structure is quite a complicated phenomenon. Ngo, Mendis, Gupta, 
and Ramsay (2007) have provided a recent overview of the problem. 


At any point in the path of an explosive blast wave, the pressure rises rapidly to a maximum value that is 


known as the overpressure. This is followed by a slower — although still quite rapid - decay of the 
pressure to below the ambient pressure. During this decay to ambient pressure, an impulse is delivered to 
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any object experiencing the increased pressure. This portion of the blast wave is then followed by a 
smaller drop in pressure below the ambient pressure and a more gradual return to atmospheric pressure. 
This portion can exert suction on an object, and can also deliver debris to the object that has been sucked 
into the blast wave. 


In a crude analysis, the effect of a blast wave on a structure is related to the magnitude of the 
overpressure. Thus, in FEMA and DOD documents, tables similar to Table 1 are often displayed. Table 
1, taken from the DOD document TB 700-2 (cited in Jeremic and Bajic, 2006), lists the types of damage 
to be expected from various blast overpressures. 


Table 1. Expected shock wave effects on objects (from DOD’s TB 700-2, cited in Jeremic and Bajic, 2006). - psi 
units added 


No. Overpressure Expected damage 
kPa psi 
1 1.0-1.5 0.15-0.22 Window glass cracks 
2 3.5-7.6 0.51-1.1 Minor damage in some buildings 
3. 7.6-12.4 1.1-1.8 Metal panels deformed 
4. 12.4-20 1.8-2.9 Concrete walls damage 
5. Over 35 Over 5.1 Wooden construction buildings demolition 
6. 27.5-48 4.0 -7.0 Major damage on steel construction objects 
ve 40-60 5.8-8.7 Heavy damage on reinforced concrete buildings 
8. 70-80 10-11.6 Probable demolition of most buildings 


However, the duration of the pressure pulse in the blast also plays a role in determining damage, i.e. the 
damage is related to the impulse as well as to the overpressure. An example of this is given in Figure 5, 
based on a figure from DOD’s TM15-1300 (cited in Ngo, Mendix, Gupta, Ramsay, 2007). Figure 5 
shows a design chart for a tempered glass panel. It shows clearly that the survival of a panel depends not 
just on the peak pressure but also on the blast impulse to which it will be subjected. 


The reason why both overpressure and impulse are important in determining structure damage becomes 
evident in a simple model that is sometimes used to describe blast wave/structure interaction: In this 
model - the so-called single degree if freedom (SDOF) model - the structure (or an element of the 
structure) is replaced by an equivalent system of one concentrated mass and one weightless spring that 
represents the resistance of the structure to deformation. 
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Peak blast pressure 
(units unspecified below) 


1 10 100 1000 
Duration of blast pressure (msec) 


Figure 5. A design chart for a tempered glass panel (TMS-1300, Figure 6-43). 
It displays the peak pressure that different panels can withstand as the blast duration is varied. 


It is well known that the response of a mass on a spring to a transient force is qualitatively different when 
the natural oscillation period (determined by the spring constant and the mass) is greater than the force 
duration than when it is less than the force duration. In the former, the impulse (the integral of the force 
over time) determines the initial velocity of response of the mass, while in the latter the displacement of 
the mass is practically proportional to the force throughout its application. In a building, different 
elements of the building will have different natural oscillation frequencies, and the response of any 
particular element in its surroundings will respond accordingly either to the blast overpressure exerted on 
it or to the blast impulse. 


In Ngo, Mendix, Gupta, Ramsay, 2007, a half dozen computer programs are listed that are used for 
calculating structural responses and for designing blast-resistant structures. 


4. Scaling laws for blast overpressure and impulse 

Several approximate scaling laws have been suggested in the literature to relate blast overpressure and 
impulse - the variables of interest for determining structural damage — to the TNT-equivalent mass and 
standoff distance of the explosive. Several examples are listed below: 

Overpressure 


¢ M.A. Sadovski [cited in Jeremic and Bajic, 2006] 


p(m,r) = 0.085(m'7/r) + 0.3(m??/r?) + 0.8(m/r?) MPa [6] 
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° HLL. Brode (1955) 
p(m,r) = 0.1 + 0.67 (m/r*?) MPa for p(m,r) > 1 MPa [7a] 


p(m,r) = 0.0975(m'7/r) + 0.1455(m7/r’) + 0.585(m/r*) — 0.0019 MPa 


for 0.01 < p(m,r) < 1 MPa [7b] 
°  C.A. Mills (1987) 
p(m,r) = 0.108(m""/r) - 0.114(m77/r’) + 1.772(m/r*) MPa [8] 
°  M. Held (1983) 
p(m,r) = 2(m77/r’) MPa [9] 


Impulse 
¢ M.A. Sadovski [cited in Jeremic and Bajic, 2006] 
i(m,r) = 200(m7"/r) Pa-s [10] 
°¢ M. Held (1983) 
i(m,r) = 300(m7"/r)_ Pa-s [11] 


Note: In these expressions, the TNT-mass equivalent m is expressed in kg and the standoff distance r is 
expressed in meters 


It is interesting to see how varied the scaling laws for overpressure are, and how not many scaling laws 
have been proposed for the impulse. This is a reflection of the complex nature of the explosive blast 
wave. 


For example, the impulse depends strongly on the shape of the decaying pressure pulse following the 
large overpressure due to the blast shock wave. Multiplication of the overpressure by cited durations for 
the pressure pulse do not give good approximations for the impulse. 

In the following we shall use Sadovski’s expressions for both the overpressure and the impulse, i.e. we 
shall use eq. [6] for p(m,r) and eq. [10] for i[m,r]. Sadovski based his scaling laws on numerous 
experimental results; and in addition, for the range of overpressures of interest, comparisons of the 
shapes of actual data curves with curve fitting expressions seem to favor expressions containing multiple 
terms [See, e.g., G.F. Kinney (1962)]. 

Equations [6] and [10] are rewritten below in units of psi, psi-msec, lbs, and ft 


m,r) = 31.11(m"?/r) + 276.9(m77/1") + 1863(m/r’) psi form inlbsandrinft [6] 
p 


i(m,r) = 56.25(m7"/r) psi msec for m in lbs and r in ft [10] 
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The blast overpressure p(m,r) given by eq. [6] as a function of TNT equivalent mass m and standoff 
distance r is displayed in Figure 6. Figure 7 shows the blast impulse i(m,r) given by eq. [10]. 


Overpressure psi 


€0000 


Figure 6. Blast overpressure p(m,r) as a function of TNT equivalent mass m and standoff distance r, from eq. [6]. 
The z-axis shows the overpressure p(m,r) in psi. The x-axis shows the TNT-equivalent mass m in lbs, and the y-axis 
shows the standoff distance in feet. 


Impulse psimsec 


Figure 7. Blast impulse i(m,r) as a function of TNT equivalent mass m and standoff distance r, from eq. [10]. The 
z-axis shows the impulse i(m,r) in psi msec. The x-axis shows the TNT-equivalent mass m in Ibs, and the y-axis 
shows the standoff distance in feet. 


It is also of interest to solve eqs. [6] and [10] for TNT-equivalent mass m in terms of blast overpressure p 
and impulse i, as well as for standoff distance r in terms of p andi. The solutions for m(p,i) and r(p,i) are 
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rather messy algebraically, but can be obtained straightforwardly with Mathematica symbolic algebra 
software. This yielded three possible roots each for m(p,i) and r(p,i) where two of the three are complex 
and therefore nonphysical. 


The dependence of the standoff distance r(p,i) on overpressure p and impulse i is shown in Figure 8, 
whereas Figure 9 shows how the TNT-equivalent mass m(p,i) depends on p and 1. 


Figure 8. Standoff distance r as a function of overpressure p and impulse i. The standoff distance is shown along 
the z-direction in feet, with the overpressure shown along the x-axis in psi and the blast impulse shown along the y- 
axis in psi msec. 


Figure 9. TNT-equivalent mass m as a function of overpressure p and impulse i. The TNT-equivalent mass is 
shown along the z-direction in Ibs, with the overpressure shown along the x-axis in psi and the blast impulse shown 
along the y-axis in psi msec. 
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The plots show that the largest values for standoff distance r and TNT-equivalent mass m occur when the 
overpressure, p, is small and the impulse, i, is large. 


5. Design curves for blast protection 


In order to obtain an estimate for the degree of blast protection, the probability considerations need to be 
combined with a design curve that sets the limits of integration for the probability function. 


Simple SDOF model For calculating the response of a structure such as a window to a blast wave, an 
approximation that is sometimes used is to replace the components of the structure by a simple mass M 
held in place by a spring of spring constant K. This approximation is termed the single degree of freedom 
(SDOF) approximation. 

The displacement y(t) of the mass in response to a time dependent force F(t) is obtained from the equation 


d’y/dt? + w’y = F/M [11] 


where the natural oscillation angular frequency o is given by w= K/M. If the force is due to a transient 
pressure pulse P, we can further write F= PA where A is the area of the component. 


Suppose we characterize the transient pressure pulse by the crude approximation 

P(t) = Pexp(-t/T) [12] 
The impulse associated with this is I = PT 
Then, assuming that the initial displacement and velocity of the component are zero, the solution is 
simply 

y(t) = afexp(-t/T) -cos(wt) +(1/(@Tsin(@t))] [13] 
where 

a. = (PA/M)I?/( P?+(ol)”) [14] 
The design curve can be obtained by setting the maximum value of y(t), i.e., where dy/dt = 0, equal to a 
critical displacement y,, sometimes taken as 1/175 of the linear dimension of the component 
[http://en.wikipedia.org/wiki/Curtain_wall]. 
Design curve equations in dimensionless variables This then leads to the equations: 

exp(-xp/i) = (1/p) +cosx [15] 

sinx = 1/i [16] 


where the dimensionless variables p and i have been introduced. They are related to the actual 
overpressure P and impulse I by 


p= Polax 17] 
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i=I/ ax [18] 
where 

X% = yeM/A [19] 
The solution of eqs. [18] and [19] are shown in Figure 10. It was obtained by numerical solution. 


Pp 


1.5 2 2.5 3 3.5 
i 
Figure 10. SDOF design curve obtained from eqs. [18] and [19] for the dimensionless variables p and i. 
The dimensionless variables are related to the actual overpressure and impulse by p = P/w’x ,i=I/ wy, 


where x, = y-M/A is related to the critical displacement, mass, and component area. 


The design curve displays two asymptotes, one for p>>1 and i approaching 1, and the other for i>>1 and 
p approaching 4. It is easy to see from eqs. [18] and [19] that the asymptotic behaviour is described by 


p= [l-awy? when p>>1 [20] 
p=%+ (w/8)(1/) when i>>1 [21] 


The design curve of Figure 10 in the dimensionless variables p and i is a universal curve: it applies 
whenever a structural component can be represented by a SDOF mass attached to a spring. 


Design curve in actual overpressure and impulse To relate the design curves in the dimensionless 


variables to the design curves for an actual situation, the natural angular oscillation frequency of the 
system and the parameter y must be known. 


As an example, suppose we take the natural period of the system to be 0.11 sec. (@ = 56 sec"). This is in 
the range of the values shown in a specific report provided by an explosives consultant. In addition, let us 


take the curve to pass through the design point: 8 psi and 170 psi msec. 


From eqs. [20]&[21] we see that for a design point (Pu,Iu), the corresponding point (pg,ig) in the 
dimensionless space satisfy the relation 


Pa/la = @ (pa/ia) [22] 
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With the choice of the natural oscillation frequency of 56 sec” along with the design point (Py= 8 psi and 
Iy = 170 psi msec = 0.17 psi sec), eq. [25] requires pg/ig = 0.84. From Figure 10, it can be seen that pa = 
1.01145 and ig = 1.2. For this example, then, 


x = Ia/ig © = 0.0025 [23] 
Accordingly 

P=8 psi [24] 

I= 142 psi-msec [25] 


This design curve is shown in Figure 11. 


P(psi) 


4U 


30 


100 150 200 250 300 350 
I(psi-msec) 


Figure 11. Design curve passing through 8 psi and 170 psi msec for a natural oscillation angular frequency of 56 
1 
sec. 


As another example, suppose instead the design point is taken to be 16 psi and 185 psi msec, and the 
natural oscillation angular frequency remains at 56 sec’. The same procedure leads in this case to 


P=9.42p psi [26] 
I=1681 __ psi-msec 


The resulting design curve for 16 psi and 185 psi msec is shown in Figure 12 (dashed) and compared to 
the design curve of Figure 10 for 8 psi and 170 psi msec. 
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Pressure (psi) 


100 150 200 290 300 350 


Blast impulse (psi msec) 


Figure 12. Design curves for blast protection. The solid curve is chosen to pass through 8 psi and 170 psi msec for a 
natural oscillation angular frequency of 56 sec’. The dashed curve is chosen to pass through 16 psi and 185 psi 
msec, with the same natural oscillation angular frequency of 56 sec’. 


Design curve translated into (r,m) space The level of blast protection depends on the particular design 
curve chosen. In the next section, the design curves will be combined with the probability function 
discussed in Section 2 to obtain a confidence level for blast protection. For that purpose, it will be 
convenient to transform the design curves depicted in Figure 12 in the (overpressure, impulse) plane to 
the (equivalent TNT mass m, standoff distance r) plane. This can be done using eqs. [12] and [13]. 
Figure 13 displays the design curves of Figure 12 in the (m,r) plane. 


Equivalent TNT 
mass m (lbs) 


190 20 20 30 330 40 453 50 


Standoff distance r (feet) 


Figure 13. Design curves for blast protection. These are the same design curves as in Figure 12, but plotted here in 
the (m,r) plane. The solid curve is chosen to pass through 8 psi and 170 psi msec for a natural oscillation angular 
frequency of 56 sec’'. The dashed curve is chosen to pass through 16 psi and 185 psi msec, with the same natural 
oscillation angular frequency of 56 sec’. 
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6. Probability of blast protection 


The significance of a design curve in Figure 13 is that the associated structural component is protected 
from any blast for which the TNT equivalent mass and standoff distance gives a point lying under the 
curve. Accordingly, the probability that the component will withstand a blast is obtained by integrating 
the probability function over all points (m,r) that lie underneath the curve. 


The design curves of Figure 13 are almost linear through the portion of the region they occupy where the 
probability function is zero. Linear fits to the curves in this portion are shown superimposed on the 
contour diagram for the probability function (of Figure 4), in Figure 14. The curves are shown only in the 
range of m and r where the probability function is nonzero. 


Equivalent TNT 
mass (Ibs) 


300 350 400 450 500 


Standoff distance (feet) 


Figure 14. Design curves of Figure 13 superimposed on the probability function contour diagram of Figure 4. 
Again, the solid design curve is chosen to pass through 8 psi and 170 psi msec for a natural oscillation angular 
frequency of 56 sec”, whereas the dashed curve is chosen to pass through 16 psi and 185 psi msec, with the same 
natural oscillation angular frequency of 56 sec’. The probability function F,,(m,r) is given by eq. [5] as the product 
of truncated normal distributions in m and r. 


For both design curves, it is seen that the maximum of the probability function lies below design curves. 
It is also seen that more of the probability function is below the dashed curve (16 psi, 185 psi msec 
design) than below the solid curve (8 psi, 170 msec). Thus, it would be expected that the confidence 
level for protection for the former would be higher than for the lower. 


For the design curves depicted in Figure 14, the probability Pprotectea that the structural component will be 
protected against a blast is 


Ima Smaller of {Mecsign(T), Mmax} 


P protected = Sdr SFme(m,r) dm [27] 


Tin Min 
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where Mdesign(t) denotes the design curve. For the two representative design curves of Figure 14, this 
integration gives: 


Pprotected (8 psi, 170 psi msec, 56 sec’ design) = 80% 
Prrotected( 16 psi, 185 psi msec, 56 sec’ design) = 91% 


i.e., the probabilities of blast protection for the two designs are 80% and 91%, respectively. 


7. Summary and discussion 
7a. Summary 


A probabilistic risk assessment procedure has been followed to estimate the level of blast protection at a 
facility. We have assumed that the damage would be collateral, with the principal terrorist target being a 
nearby facility. A probabilistic approach was adopted because there are uncertainties associated with 
both the TNT-equivalent mass of the explosive and with the standoff distance of the detonation. These 
uncertainties in explosive mass and standoff distance translate through commonly used scaling laws into 
uncertainties in the blast wave overpressures and impulses. The overpressures and impulses are of more 
direct interest to structural designers. 


Normal distributions, with cutoffs at reasonable minimum and maximum values, have been assumed for 
both the TNT equivalencies of the charge and the standoff distances [Figure 3]. A structural component 
design curve for blast mitigation has been derived from an oft-used simple single degree of freedom 
(SDOF) model. This design curve defines the boundary between the overpressures and impulse 
combinations from which the structural component will be protected and those for which the design 
affords no protection. Use of the scaling laws translates this information into a boundary between 
combinations of TNT-equivalent mass and standoff distance. The probability of protection is then 
obtained by integrating the normal probability distributions over all masses and distances that lie on the 
protected side of the design curve. 


Two numerical examples have been given, both for a structural component that has a natural oscillation 
period of 0.11 sec, corresponding to a natural oscillation angular frequency of 

56 sec’. In the first example, the design curve has been chosen to pass through the point (8 psi, 170 psi 
msec); the resulting probability of protection is 80%. In the second example, the design curve has been 
chosen to pass through the point (16 psi, 185 psi msec): the resulting probability of protection for this 
case is 91%. Thus, the first example gives protection at the 1.3-sigma level, whereas the second example 
provides protection at the 1.7-sigma level. 


7b. Discussion 


This report has focused on deriving the probability of protection afforded by a specific window design as 
a result of an explosive blast. 


It should be noted that as more information is known about risk mitigation such as the emplacement of 
bollards to enforce standoff, the geometry of a target facility relative to the surrounding streets, and likely 
screening procedures, improvements can be made on the normal truncated probability distributions 
assumed for TNT-equivalent masses and standoff distances. 
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In addition, we believe that the probabilistic approach in this report and applied to a specific terrorist- 
initiated explosive threat can find application in a much broader range of operational risk problems. A 
general four-step process is envisioned: 


1. The starting point is to assume an appropriate probability function (often a normal distribution) for a 
set of parameters that is inherently random, and that have some (indirect) connection with the mitigation. 


2. The second step is to develop equations (scaling laws) that relate these inherently random parameters 
to a second set of parameters that most directly impact the mitigation under consideration. 


3. The third step is to identify design curves (or surfaces) for a group of selected risk-aversion measures 
in the second set of parameters. 


4. The fourth step is to use identified scaling laws to map the design curves (surfaces) onto the space of 
inherently random parameters, and to integrate the normal probability function over that portion of the 
space that is protected by the mitigation design. This then gives the probability that the applied 
mitigation will be successful, i.e. it provides statistical confidence in a successful implementation. 


Alternatively, steps 3 and 4 can be replaced by: 


3 (alternate). The third step is to use the results of the first two steps to generate a probability distribution 
for the second set of parameters. 


4 (alternate). The fourth step is to integrate the probability function over the portion of the space of the 
second set of parameters that is defined by the design curves (surfaces) for selected risk mitigation 


measures. This then gives the probability that the risk mitigation will be successful. 


The choice between Steps 3 and 4 or the alternate Steps 3 and 4 could be based on the relative ease of 
calculation for the two options in a given risk scenario. 


The approach can be used to systematically estimate the effectiveness of mitigation measures and 
associated costs required to achieve varying degrees of protection. 
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